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ABSTRACT 

We show that interferometry can be applied to study irregular, rapidly rotating structures, as are 
expected in the turbulent accretion flow near a black hole. Specifically, we analyze the lagged covari¬ 
ance between interferometric baselines of similar lengths but slightly different orientations. For a flow 
viewed close to face-on, we demonstrate that the peak in the lagged covariance indicates the direction 
and angular velocity of the emission pattern from the flow. Even for moderately inclined flows, the 
covariance robustly estimates the flow direction, although the estimated angular velocity can be sig¬ 
nificantly biased. Importantly, measuring the direction of the flow as clockwise or counterclockwise 
on the sky breaks a degeneracy in accretion disk inclinations when analyzing time-averaged images 
alone. We explore the potential efficacy of our technique using three-dimensional, general relativistic 
magnetohydrodynamic (GRMHD) simulations, and we highlight several baseline pairs for the Event 
Horizon Telescope (EHT) that are well-suited to this application. These results indicate that the EHT 
may be capable of estimating the direction and angular velocity of the emitting material near Sgr A*, 
and they suggest that a rotating flow may even be utilized to improve imaging capabilities. 

Keywords: accretion, accretion disks - black hole physics - Galaxy: center - techniques: high angular 
resolution - techniques: interferometric 


1. INTRODUCTION 

The innermost accretion flows around black holes are 
the subject of intense numerical study despite a dearth 
of observational constraints. Even the Galactic Center 
supermassive black hole, Sagittarius A* (Sgr A*), has 
yielded only limited concl usions about its accretion en- 
vironment and dynamics (Genzel et al. 2010[ Yuan & 


Narayan]|20I4 ), with x-ray observations constraining the 


accretion boundary conditions on scales comparable to 
the Bondi radius, at roughly 10® times the gravitationa l 
radius re = GM/c^ of Sgr A* (Baganoff et al. 2003), 


and with radio observations providing rough estiniates of 


ro) at wavelengths from 3 — 13 mm ( 

Lo et al.|1998 Shen 

et al. 120051 IBower et al.|2006| Lu et aL] 

wrn 

IBower et al. 

2014 Cwinn et al. 2014|). Moreover, 

because the radio 


tering, it was measurements of polarization and Faraday 
rotation that conclusivel y constrained the flow properties 


on these smaller scales (Aitken et al. |2000 

Quataert & 

Gruzinov||2000 

Bower et al.|2003 

Marrone et al.|2007|. 

This limited 
abruptly with t 
scope (EHT), e 
0.87-mm very-- 

observational perspective will change 
he completion of the Event Horizon Tele¬ 
project to develop a global 1.3-mm and 
ong-baseline interferometry (VLBI) net- 


work (Doeleman et al. 2009a). Ultimately, this net¬ 
work will provide a nominal angular resolution of tens 
of microarcseconds, sufficient to resolve the event hori¬ 
zons of the nearest supermassive black holes, including 
the Galactic Center black hole, Sgr A*. The EHT can 
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achieve detections on Sgr A* with an integration time of 
minutes or less - significantly shorter than the period of 
the innermost stable circular orbit of Sgr A* - rendering 
it sensitive to both steady features and variability in the 
emission from Sgr A*. 


Past EHT observat ions of Sgr A* (Doeleman et al. 
2008 Fish et aL]|201I ) have already suggested structure 


that is more compact than the size of the photon ring 


that bounds the black hole “sha dow” (Bardeen 1973 Fal- 
cke et al.|2000 Takahashi|2004 ). This compact structure 


IS most commonly explained via an accretion disk with 
its angular momentum axis inclined relative to the line 
of sight, which results in an image that is dominated by 


the disk (e.g 

.jlDexter et al. 

2010 

Broderick et al. 

2011 

Psaltis et al. 

2(J151. However, given that the current data 

are extreme. 

y sparse, the inclination and image proper- 


ties cannot yet be confidently constrained. In addition, 
current estimates of the inclination are subject to a de¬ 
generacy between supplementary inclinations, {9, n — 0}, 
(see Figure because simulated images of accretion 


tion axis (IMoscibrodzka et al. 

2009 Dexter et al. 

2010 

Broderick et al. 1120111 |Shcher 

oakov et al.||2012| |P 

saltisi 

et al.||2015D. 


it is equally important to develop model-independent as¬ 
sessments of the data. Chief among these is syn thesis 
imaging (see, e.g., Lu et al.||2014[ Fish et al.|| 2014|). Yet 


2008 


Gillessen 


with a mas s of 4 x 10® Mq (Chez et al. 
et al.|200'9 ), Sgr A* has a gravitational timescale of only 
GM/c® « 20 seconds and an orbital period at the inner¬ 
most stable circular orbit of only 4—30 minutes, depend¬ 
ing on spin (Bardeen et al.|1972 ). These short timescales 
suggest that conventional Earth-rotation synthesis imag¬ 
ing will be inapplicable for Sgr A* and, in the best case. 
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will ignore the rich physics encoded within the variability, 
such as turbulence, orbital motion, and flaring. Thus, in 
contrast with past work to infer properties of the quies¬ 
cent image using non-imaging EHT data products, our 
present emphasis is to study dynamics of the emitting 
material with EHT data. 

Several authors have already explored how non¬ 
imaging VLBI can be a pplied to study rapid temporal 
variability. For instance Broderick & Loeb (2005, 20061 


simulated an orbiting “hot spot” around Sgr A* and cal- 
culated the expected images at submil l imeter and near- 


infrared wav elengths. Doeleman et al. (2009b) and Fish 


et al. (20091 then showed that EH'i' baselines can sensi¬ 


tively detect periodicities associated with these hot spots. 

bi 


Addre ssing more general circumstances, [Broderick et al. 
(|201l|) suggested that phase-referenced observations with 
the EHT may allow microarcsecond tracking of the im- 
age ce ntroid on orbital timescales, and [Johnson et al.j 
(2014) showed that polarimetric VLBI with the EHT is 
capable of microarcsecond astrometry of compact flaring 
structures, even for faint, non-periodic flares. 

In this paper, we explore a different metric: the tempo¬ 
ral covariance between pairs of baselines. For baselines 
that are of a similar length and orientation, this covari¬ 
ance is a sensitive probe of image rotation. Because real¬ 
istic flows around black holes are subject to shearing and 
other secular evolution, these measurements can be used 
to determine orbital timescales using nearby “snapshots” 
between which the flow undergoes stable rotation. The 
temporal covariance then reflects the direction of the flow 
on the sky (clockwise or counter-clockwise), the angular 
velocity of the emitting material, and the radial distri¬ 
bution of emitting material. Estimating the direction 
of the flow immediately breaks the degeneracy in sup¬ 
plementary inclinations when analyzing quiescent images 
permitting an unambiguous determination of the angular 
momentum axis of the accretion flow. Importantly, our 
work requires no assumptions about the appearance of 
the quiescent structure or about precise flow dynamics. 
We only assume that there is an inhomogeneous, rotating 
component. 

We begin, in with a brief discussion of accretion and 
orbital dynamics near a black hole. Next, in ® we derive 
interferometric relationships for a differentially rotating 
flow. Then, in 0 we apply our technique to synthetic 
observations of three-dimensional GRMHD simulations. 
In ^ we consider practical limitations when using this 
technique on EHT data products and discuss the EHT 
baseline pairs that are especially well-suited to this pur¬ 
pose. Finally, in we summarize our main results. 


2. THE ACCRETION FLOW NEAR A BLACK 
HOLE 

The perceived angular velocity of a rotating accretion 
flow by a distant observer depends on properties of the 
spacetime, the accretion flow itself, and the viewing ge¬ 
ometry. We now discuss each of these in turn. 


2.1. Effects of the Spacetime on Orbital Velocities 

In the Kerr spacetime for a rotating black hole with 
dimensionless spin 0 < a < 1 and gravitational radius 
re = GM/c^, the period Porb of a circular equatorial 
orbit at radius r (in Boyer-Lindquist coordinates) can be 



Figure 1. Illustration of the degeneracy in supplementary in¬ 
clinations of an accretion disk. For accretion disks with the two 
inclinations shown, the Doppler boosting will be identical (neglect¬ 
ing effects from black hole spin). As a result, the two images will 
have similar time-averaged appearances, but the direction of the 
flow on the sky will be opposite. An extreme example of this de¬ 
generacy is an accretion disk viewed face-on. 


written as ( Bardeen et al.|[l972 ) 

PoTh{r) = 27r {r/reff/'^ ± a 




( 1 ) 


where +/— corresponds to prograde/retrograde orbits. 
For a Schwarzschild black hole (a = 0), this expression 
reduces to Porb(^) = 2tt s/r^JjGM), which is familiar 
as Kepler’s third law. Because of the differential rota¬ 
tion implied by Eq. large coherent features in a Keple- 
rian flow will shear apart on an orbital timescale, quickly 
eliminating pure periodicities in light curves]^ For flows 
viewed off-axis, the effects of strong gravitational leasing 
can also change the apparent angular velocity depending 
on orbital phase. 

Although no significant periodicities have been de¬ 
tected for Sgr A*, some simulations suggest that the 
radio emission is dominat ed by material at a Boyer- 
Lindquist radius of ~5rG (Shiokawa 2013). Because all 
emissions from small radii are leased to similar apparent 
radii (see, e.g., Broderick et al.|2009l, time-averaged im¬ 
ages may not be able to constrain the emission radius as 
tightly as the orbital dynamics. However, orbital period¬ 
icities will probably be unable to meaningfully constrain 


the b lack hole spin, as noted by Broderick & Loeb (2005 
2006). For instance, at a radius of dre, the effect of spin 


on the period of prograde orbits is <10% and is degener¬ 
ate with a 6% change in emission radius. Thus, given the 
uncertainty of the precise emission radius at submillime¬ 
ter wavelengths, the black-hole spin can probably not be 
securely estimated from its effect on orbital periodicities. 

However, while the orbital period at a given radius is 
rather insensitive to spin, the radius risco of the inner¬ 
most stable circular orbit (ISCO) is not. For a maxi¬ 
mally spinning black hole, risco ranges from cq to Ore 
for prograde and retrograde orbits, re spectively. This 
steep dependence on spin was used by [Doeleman et al.j 
(20121 to infer that the accretion disk in M87 is undergo- 


ing prograde rotation by associating the apparent size of 
the emission region with the lensed ISCO. Because mea¬ 
surements of the orbital period can accurately estimate 
the emission radius, these estimates could then provide 


For instance, in a purely Keplerian flow Porb (r(l + /)) ~ 
[!+§/] Porb (^) 5 SO a feature situated at a radius r and having 
a radial extent of 0.16r will shear by approximately a quarter of 
an orbit each orbital period. 
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a meaningful upper-bound on risco and, thus, a lower- 
bound on the signed spin (-1-/— for prograde/retrograde). 

2.2. Effects of the Accretion Flow on Orbital Velocities 

Because Sgr A* is highly under-luminous, emitting 
only of its Eddington luminosity, the accretion 


flow near Sgr A 

* is most likely a hot, thick disk (e.g.. 

Yuan et al. 120031 

Loeb & Waxman|2007| 

|Yuan & Narayan] 

2014). Such disks tend to have suh 

tion profiles because of their strong 

)-Keplerian rota- 
pressure support 

(Narayan & Yi||1994 

Yuan & Narayan 

2014). The frac- 


rian fix is second order in the fractional scale he ight H/R 


of the disk (e.g., Blandford & Begelman 1999 Eq. 20). 
GRMHD simulations ot bgr A*" typically have H /R rang¬ 


ing from 0.1 to 0.3 (e.g., Dexter et al. 2010), so the 
pressure support may decrease the angular velocities by 
<10%. For the particular GRMHD simulation used in 
this paper (see (Q, the average azimuthal fluid velocity 
differs from the Keplerian velocity by only ~1% outside 
the ISCO. 

A potentially more serious limitation concerns the pat¬ 
tern velocities of emission features, which may differ from 
their underlying fluid velocities^ For the GRMHD sim¬ 
ulations that we discuss in ® the pattern velocity is 30- 
40% lower than the fluid velocity at the radius of peak 
emission (wSro); the difference becomes even more pro¬ 
nounced closer to the event horizon. Although we are 
not aware of simulations that exhibit a predominantly 
counter-rotating submillimeter emission pattern, the re¬ 
lationship between the pattern and fluid velocities is an 
important consideration when interpreting measured an¬ 
gular velocities and will be analyzed in depth elsewhere 
(Shiokawa et ah, in preparation). 


2.3. Effects of the Viewing Geometry on Orbital 
Velocities 


When rotation is not viewed face-on, apparent angular 
velocities depend on orbital phase. For instance, consider 
a point particle in a circular orbit at an inclination 0 < 
0 < 7r/2 relative to the line of sight (i.e., the motion of 
a bead on a circular wire). At an orbital phase (j), the 
apparent angular velocity is 


9) oc 


COS0 

cos^ 6 + cos^ (j) sin^ 9 


( 2 ) 


At ((i = 0 the particle is moving away from the observer, 
and at 0 = TT is moving toward the observer. The veloc¬ 
ity component along the line of sight is cos(()sin0. Note 
that the velocity is zero for an edge-on view of the orbit 
{9 = 7r/2) because in that case the motion is restricted 
to a line and has no angular velocity. Near a black hole, 
relativistic aberration and leasing will also affect the ap¬ 
parent velocity. For instance, emission will be lensed 
above and below the black hole, producing a non-zero 
apparent angular velocity even for the edge-on case. 

The orbit-averaged angular velocity must, of course, 
agree with the true angular velocity. However, because 
of Doppler effects and leasing, the emitting material will 
vary in brightness through the orbit. For instance, let¬ 
ting V = 7“^(1 —/3||)“^ denote the Doppler factor, where 


® We thank the referee for identifying the importance of the 
emission pattern velocity. 


7 = (1 — is the Lorentz factor and /3 is the nor¬ 

malized velocity in the rest frame of the observer, the 
observed flux density is scaled by a factor of 
ative to that in the co-moving fra me, where a is the 
spect ral index of the emission (e.g., Blandford &: Konigl] 


1979). Gonsequently, the orbit-averaged angular velocity 


weighted by VA will underestimate the true angular ve¬ 
locity; most emission comes when the orbit is approach¬ 
ing the observer, with a low apparent angular velocity. 

An additional complication arises because the angular 
velocity must be defined relative to a particular centroid 
- typically the centroid of the quiescent flux. For an 
inclined flow, because of Doppler boosting on the ap¬ 
proaching side of an accretion disk, this centroid will not 
be centered on the black hole. We will now derive a 
strategy for estimating the angular velocity using inter¬ 
ferometric visibilities that mitigates this effect. 


3. LAGGED INTERFEROMETRIG COVARIANGE 
3.1. Interferometric Observations of a Rotating Flow 


We now explore how signatures of a rotating flow are 
manifest in interferometric observables. The interfero¬ 
metric visibility /(u) measured by a baseline u is re¬ 
lated to the source brightness distribution /(x) via the 
van Gittert-Zernike Theorem (Thompson et al. 2001): 


i(u,t)= / cfixl{x,t)e 


( 3 ) 


In this expression, u is the vector baseline orthogonal 
to the line of sight, in wavelengths, and x is an angular 
coordinate on the sky, in radians. We have included a 
time coordinate, t, to account for the possibility of a 
changing source emission structure with time. 

From Eq. it is evident that an image rotation by 
some angle 9 about the origin x = 0 leads to an identical 
rotation in the visibility domain about u = 0. Moreover, 
standard interferometric observables - visibility ampli¬ 
tudes and closure phases - are unaffected by a shift of 
the image center: x —>• x-f Axg. Gonsequently, these ob¬ 
servables for one set of baselines sampling the unrotated 
image will be equivalent for the same baselines rotated 
by 9 but sampling the image after it is rotated by 9 about 
any fixed point Xq. This property helps to mitigat e the 
angular-velocity bias from an offset image center ( p.3| ) 
when studying angular velocities. 

Hence, as long as the image is not azimuthally sym¬ 
metric, stable rotating structures will introduce a lagged 
correlation between visibilities on pairs of baselines with 
similar lengths but different orientations. The angular 
velocity of the rotating flow is then given by the angu¬ 
lar difference between the baseline directions divided by 
the temporal lag corresponding to the peak covariance. 
This inference, which determines both the direction and 
angular velocity of rotation, is determined entirely in the 
visibility domain and can be achieved with as few as two 
baselines (three stations). 


3.2. Identifying the Peak Lagged Correlation 

We now describe our procedure to estimate the lagged 
correlation. To properly identify the peak lag, we 
must address three potential sources of contamination: 
1.) There will be overall changes in the total flux of the 
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Time (Minute) 


Figure 2. Left panels show normalized visibility as a function of time for a simulated GRMHD movie sampled on three baselines of 
identical lengths (8 GA) but slightly different orientations. Results are shown for four viewing inclinations: 9 = 0 ° (face-on), 9 = 30°, 
9 = 60°, and 9 = 85°. Corresponding images on the right show single frames for e ach case, with a linear color scale. For these images 
and visibilities, we show results before accounting for interstellar scattering (see ^5.4[ |. At low inclinations, the lags between the time series 
are readily apparent by comparing the turning points of each curve and arise because of rotation of the irregular flow. Especially for the 
face-on disk, these lags accurately estimate the orbital periodicity, even though the visibility curves show no clear periodicities. In contrast, 
for the nearly edge-on case (9 = 85°), the lags are irregular in both magnitude and sign. 


image (i.e., the zero-baseline visibility), 2.) There will 
be slow secular evolution of the bulk emission structure 
and of observing parameters (e.g., from rotation of the 
Earth), and 3.) The brightness centroid may not be cen¬ 
tered on the black hole. 

The first consideration is especially important for short 
baselines, where variations in the interferometric visibil¬ 
ity will be tightly correlated with modulation of the total 
image flux. As a result, lagged correlations will have a 
peak at zero lag. To eliminate this feature and to ac¬ 
count for source flux modulation, it is advantageous to 
work with normalized visibilities - i.e., visibilities divided 
by the simultaneous zero-baseline visibility. To mitigate 
the second contaminating effect, we subdivide each long 
time-series into shorter segments and determine the peak 
lag for each separately. This subdivision also naturally 
accommodates observational constraints, such as regu¬ 
lar breaks for calibration or pointing scans. Finally, we 
eliminate the third source of contamination by studying 
only the normalized visibility magnitudes, as discussed 



In each segment, we then estimate the lagged cross¬ 
correlation p{At) using the classical estimator: 


p{At) = 


{[A{t) - {Am [B{t+At) - {B{t+At))]) 


(4) 


where {A{t),B{t)} are the normalized visibilities 
/(u, t)//(0, t) on the pair of baselines, and ax denotes 

the standard deviation of the time series x. 

When the EHT is complete and begins collecting regu¬ 
lar data on Sgr A*, additional knowledge of the variabil¬ 
ity can be applied to develop more sophisticated estima¬ 
tors of the lagged covariance. For example, differences 
of nearby measurements (approximating the time-series’ 
derivatives) are effect ive for de-trending and whitening 
stoch astic time series (|Brockwell & Davis|2002 Box et al. 
2008|), and may facilitate superior estimates of the lag. 
Also, alternati ve metrics such as the Discrete Correla¬ 
tion Function ( Edelson fc Krolik[|1988 1 could be adopted 
for unevenly sampled data, such as from irregularly in¬ 
terspersed scans on calibration targets, although VLBI 
scans can be correlated at arbitrarily short segmenta¬ 
tion times. However, since our focus is merely a proof- 
of-concept, we will use the simple correlation estimate 
defined by Eq. 


4. APPLICATION TO GRMHD SIMULATIONS 

We tested our new technique by applying it to a 3D 
GRMHD simulation (bO-high from Shiokawa 2013) of 
a rad iatively inefficient accretion flow (e.g., Fsm et al. 
1996) onto a massive (4.5x10® Mq) spinning black hole 
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( Gammie et al.|[2003 Noble et'n^|2006 ). We started the 
simul ation with a hot, geometric ally thick, and tenuous 
disk (Fishbone & Moncrief||1976 1 around the black hole. 
The disk was seeded by a weak poloidal magnetic field so 
that the magnetorotational instability (MRI) could grow, 
driving accretion. The radii of the initial inner edge and 
pressure maximum of the disk were 12rQ and 24rG, re¬ 
spectively. We set the dimensionless spin parameter of 
the black hole to be a = 0.9375, following the “best - 
bet model” for Sgr A* from Moscibrodzka et al. (2009). 
For this spin, the ISCO radius is risco = 2.04^^ with 
a corresponding orbital period of Fisco = 24.25tG = 
8.96 minutes, where to = GM/c^ = 22.17 seconds. How¬ 
ever, most of the 230 GHz emission originates from ma¬ 
terial at a Boyer-Lindquist radius of 6 — 5rG for viewing 
inclination of 0° — 80°, respectively, wi th a correspond - 
ing orbital period of Porb 30 minutes (|Shiokawa|2013|) . 
Note that spin has a <10% effect on orbital period at 
this radius and amounts to at most a factor of 2 even at 
r = re- 

Our 260 X192 X128 simulation grid was defined by mod¬ 
ified spherical coordinates: logarithmically scaled radial 
coordinates spanning 1.22 — 240rG, poloidal coordinates 
with 2° cutouts at the poles to avoid the coordinate sin¬ 
gularities, and azimuthal coordinates spanning the full 
27r. The MRI saturated at the initial pressure maximum 
radius around t ^ SOOOtG; we then ran the simulation 
for an additional 6500tG, which defined the data used in 
our subsequent analysis. 

For the radiative transfer, we performed general rela¬ 
tivistic ray-tracing by integrating synchrotron emission 
and absorption along photons’ geodesics until they es¬ 
cape the s imulation box an d fall into each pixel of the 
“camera” (Noble et al.|2007). Unlike most previous stud¬ 
ies, our radiative transter does not use the “fast light” 
approximation, which assumes that all the photons emit¬ 
ted in one time slice arrive simultaneously at the camera. 
We instead account for evolution of the fluid as each pho¬ 
ton propagates. We assumed the electron distribution 
function to be thermal and the proton-to-electron tem¬ 
perature ratio to be 3. At A=1.3 mm, contributions from 
bremsstrahlung and Compton scattering are negligible. 
Because the disk mass can be chosen arbitrarily in the 
conversion from simulation units to physical units (the 
disk evolution is independent of its mass in the regime 
where self-gravity and radiative effects are negligible), we 
chose a value so that the simulation’s time-averaged flux 
density was com parable to observed values for Sgr A* at 


A=1.3 mm (e.g., [Bower et al. 


2015). 


To test our proposed method, we used 1700 frames 
spaced by 0.5tG ~ H seconds, equivalent to a 5.2-hour 
observation. Figure shows example image snapshots 
and time series for interferometric visibilities at three 
different viewing inclinations. Note that these time se¬ 
ries do n ot reflect the orbital periodicities, as was also 
noted by Dolence et al. (2012) for total-flux light curves 
at A=1.3 mm from similar simulations. 

To account for slow trends in the data, we divided each 
time series into 200-frame segments (37 minutes) and av¬ 
eraged the peak lags calculated separately in each, as dis¬ 
cussed in j ]3.2[ Figurej^shows the resulting peak lag as a 
function of angular baseline separation for three baseline 
lengths and for three viewing inclinations. When the ac¬ 
cretion flow is viewed face-on, the inferred orbital period 


is close to the value for material at the radius of maxi¬ 
mum emission (Porb ~ 30 minutes). At an inclination of 
9 = 30°, the inferred orbital period is only accurate for 
long baselines (>3 GA). At 9 = 60°, the inferred orbital 
periods are lower than the true value by a factor of ~2, 
even on long baselines, but the inferred direction of the 
flow on the sky is correct in every case. At 9 = 85°, the 
peak lag varies erratically in sign and magnitude, as is 
expected from the near symmetry of the image in this 
case. 

Note that when the inferred periodicity is stable (i.e., 
at low inclinations) it is faster than the orbital period 
at the ISCO for a non-rotating black hole (in this exam¬ 
ple, 34 minutes), so these measurements could be used 
to argue that the black hole spin was non-zero and that 
the accretion flow was in a prograde orbit. In practice, 
to ensure a robust measurement of angular velocity and 
direction will require confirmation at different observ¬ 
ing epochs and on different baseline pairs, ideally sam¬ 
pling different position angles. Variations among differ¬ 
ent baseline pairs and different observing epochs would 
indicate that the inferred angular velocity and direction 
are not meaningful and could provide evidence that the 
flow is being viewed at high inclination. 

5. CONSIDERATIONS FOR THE EHT 
5.1. Ideal EHT Baselines 

The EHT has several promising baselines to study 
lagged correlation for Sgr A* (see Figure |^; baselines 
from pairs of sites at similar latitude to the South Pole 
Telescope are especially well-suited to our method. For 
instance, the pair of baselines from the South Pole Tele¬ 
scope (SPT) to the Submillimeter Array (SMA) and to 
the Large Millimeter Telescope (LMT) have a mutual vis¬ 
ibility of approximately 4 hours, assuming a 15° elevation 
pointing limit. Over this entire span, the ^8 GA pro¬ 
jected baseline lengths differ by no more than 1%, while 
the angular difference between the baselines ranges from 
29.4° to 32.8°. A 30-minute orbital periodicity would 
have a corresponding peak lag of minute between 
these baselines. 

Even the current EHT may be suitable for this method. 
The baselines from the SMA to SMT and CARMA have 
lengths (~3 GA) with ratios between 0.8 — 0.9 for their 
4 hours of mutual visibility and a baseline rotation of 
5.8 — 9.2° Table [^provides details for four of the most 
promising oaseline pairs. 

5.2. Temporal Resolution of the EHT 

Because VLBI data can be correlated and sampled on 
arbitrary segmentation times, temporal resolution is not 
likely to be a fundamental limitation for our proposed 
technique. The EHT will even achieve signal-to-noise 
ratio on scans lasting only seconds or less. For exam¬ 
ple, consider the SPT-SMA baseline. This baseline will 
have an effective system equivalent flux density (SEFD) 
of ~5000 Jy. Then, for 4 GHz of bandwidth, a 10-second 
scan would have a thermal noise of ~ 20 mJy. In the sim¬ 
ulations discussed in Q the source flux ranges from 10 

® The CARMA observatory was shut down following the 2015 
EHT campaign. However, the CARMA site is still relevant for 
continued observations with the EHT because a nearby site may 
be added (associated with the Owens Valley Radio Observatory). 
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0 = 0 ° 0 = 30 ° 0 = 60 ° 0 = 85 ° 



Angular Baseline Separation 

Figure 3. Peak lag as a function of baseline separation for four inclinations: 6 = 0° (face-on), 6 = 30°, 9 = 60°, and 9 = 85° (see 
Figurej^for characteristic images). In each case, results are shown for three baseline lengths, 1 GA, 3 GA, and 8 GA (respective resolutions 
of Ri 2(KT /ras, 70 /ras, and 25 /ras). The panels also indicate the angular-velocity curves for material in circular orbits at the ISCO and for 
a 30-minute orbital period, which is close to the period for material that dominates the emission (r ~ 5rQ). For inclinations up to 60°, 
the derived direction of the flow is correct, but long baselines are increasingly important at higher inclination to measure a meaningful 
angular velocity, and an upward bias in the inferred rotational velocities is important for 6 )> 60° because of the inclination dependence 
of projected angular velocity (see, e.g., Eq. [^. For an inclination of 85° (close to edge-on), our method cannot identify a predominant 
rotation direction (by symmetry), so it produces irregular results that will vary among different baseline pairs and observing epochs. 


Table 1 

Exemplar Baseline Pairs for the EHT. 


Baseline Pair 

Mutual Visibility (hours) 

Resolution (/i.as) 

Max Length Ratio 

Angular Difference 

Lagx 30 (seconds) 

SPT-SMA/LMT 

4.0 

25 

1.01 

29.4° - 32.8° 

74-82 

SPT-SMT/GARMA 

5.0 

25 

1.02 

2.2° - 4.4° 

6-11 

SPT-PV/PdB 

2.2 

25 

1.03 

3.2° - 5.0° 

8-12 

SMA-GARMA/SMT 

3.9 

65 

1.22 

5.8° - 9.2° 

14- 23 


Note. — SPT: South Pole Telescope, SMA: Submillimeter Array, CARMA: Combined Array for Research in Millimeter-wave Astronomy, 
SMT: Submillimeter Telescope, LMT: Large Millimeter Telescope, ALMA: Atacama Large Millimeter/submillimeter Array, PV: Institut 
de Radioastronomie Millimetrique (IRAM) telescope on Pico Veleta, PdB: IRAM Plateau de Bure Interferometer. 


to 200 mJy on this baseline. Even after accounting for 
interstellar scattering, the se values will only be reduced 
by a factor of ^2 (see (5.4). Thus, the signal-to-noise ra¬ 
tio may be g, r eater than unity even on 10-second scans. 


See Lu et al. 


(2014) for a list of current EHT SEFDs. 


5.3. Effects from Calibration Uncertainties 

Precise calibration poses a major challenge for interfer¬ 
ometry and may complicate our proposed method. Al¬ 
though our method does not require phase information, 
amplitude calibration must still be stable to within the 
variability amplitude to avoid contaminating the covari¬ 
ance. This limitation may render long baselines most 
useful, where fractional variations of the signal are likely 
larger. 

Alternatively, one can eliminate stati on-based calibra¬ 
tion errors using “closure” quantities (Thompson et al. 
2001). Because the most common closure quantities 
- closure phase and closure amplitude - involve mul¬ 
tiple baselines, they are not directly useful for our 
purposes. However, fractional polarization is baseline- 
based and provides identical immunity after calibrating 
slowly changing diff erential gain and leakage terms (e.g., 
Roberts et al.|1994|) . Comparisons wi th simulations that 
include polarization information (e.g., [Shcherbakov et alT 


Shcherbakov & McKinney||2013f Uexterj|2014) will 

us to assess whether lagged correlation of fractional 


2012 

allow us to assess whether lagged 

polarization on close baseline pairs can likewise reflect 


dynamics of the bulk flow. However, because the polar¬ 
ization direction can change throughout an orbit, from 
the changing local magnetic field direction or relativistic 
aberration or from strong-held relativistic effects such as 
lensing or parallel transport, the application to polariza¬ 
tion may require signihcant modihcation. 

5.4. Effects from Interstellar Scattering 

Propagation through the turbulent interstellar medium 
scatters radio waves and causes wavelength-dependent 
blurring of images. The scattering is especially strong 
along the line of sight to the Galactic Center, and al¬ 
though interstellar scattering is subdominant to the in¬ 
trinsic structure at A = 1.3 mm, it remains an impor - 
tant consideration for the EHT (see Fish et al. 2014). 
Nevertheless, interstellar scattering is not likely to be an 
important consideration for our proposed methods. 

For instance, the dominant effect of scattering, blurring 
via a deterministic image convolution, is invertible and 
simply decreases the signal-to-noise ratio on long EHT 
baselines by a factor of up to ~4. Because the scattering 
is weaker in the North-South direction, long baselines to 
the SPT are only attenuated by a factor of ~2. Regard¬ 
less of baseline, this ensemble-average scattering effect 
will not bias the lagged covariance. 

In contrast, the subdominant effect of scattering, “re¬ 
fractive noise,” is n ot deterministic and will affect long- 
baseline properties (Narayan & Goodman 1989 Good- 
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CARMA 


SMT 


LMT 


ALMA 


SPT 



Figure 4. (Top) The planned EHT array as seen from the dec¬ 
lination of Sgr A*; PV and PdB are not visible. (Bottom) Corre¬ 
sponding baseline tracks for all planned EHT sites. Colored tracks 
show the candidate baseline pairs listed in Table^ Baseline lengths 
are given in gigawavelengths at A = 1.3 mm. Note that because 
the latitude of PV (37.1°) is nearly identical to CARMA (37.3°), 
their baselines to the SPT are nearly identical. 


man & Narayan 1989 Johnson & Gwinn 2015). However, 
retractive noise is persistent (changing over a timescale 
of >1 day) and wideband, so also will not affect the co- 
variance or dynamical imaging that we propose, which 
rely on variability timescales of minutes. 


6. SUMMARY 

We have showed that the covariance between pairs of 
interferometric baselines with similar lengths and close 
angular separation can sensitively probe the angular ve¬ 
locity of emission for a rotating flow. This non-imaging 
technique can estimate both the direction and angular 
velocity of the flow on the sky with as few as two base¬ 
lines (three stations). By employing baselines of close 
angular separation, one can accurately estimate orbital 
periods even if the rotating structures evolve significantly 
over a single orbit, as is expected from differential rota¬ 
tion in a Keplerian flow. 

Our primary motivation has been EHT observations 
of Sgr A*. While our proposed technique would be most 
effective for a face-on viewing geometry, which is disfa¬ 
vored by current VLBI constraints, we show that EHT 
baseline pairs can robustly estimate the direction and 


can roughly estimate the angular velocity even for a 
moderately-inclined rotation axis. In particular, a mea¬ 
surement of the orbital direction would break the degen¬ 
eracy in supplementary inclinations, allowing unambigu¬ 
ous comparison with larger-scale features s uch as the cir- 
cumn uclear disk, the inner stellar disk (e.g., Bartko et al.j 


2009|, and the more recently discovered Gl and G2 gas 
clouds o n trajectories passing within the Bond i radius of 
Sgr A* ([G illesseii et al.|20i2 Pfuhl et al.|20I5 McGourt] 


& Madigan 


l||2015|). 

he generali 


Despite the generality of our approach, there are signif¬ 
icant remaining uncertainties that can affect the interpre¬ 
tation of inferred angular velocities. At high inclinations, 
the estimated an^lar velocities can be significantly bi¬ 
ased (see Figure [^, and so the applicability to Sgr A*, 
which does not have a firmly established inclination, is 
not yet secure. Indeed, there is not yet a consensus on if 
the emission from Sgr A* arises in an accretion disk or a 
jet. Even for a disk viewed face-on, there may be signifi¬ 
cant differences between the pattern velocity of emission 
features and the fluid velocity. For a thick accretion disk, 
strong pressure support may affect rotational periods as 
well, invalidating direct comparisons with rotation curves 
in the Kerr metric. To resolve these remaining questions 
will require both observational input, which is imminent 
with the addition of many new EHT sites, and improved 
understanding of the relationship between accretion and 
emission properties in GRMHD simulations. 

Although we have focused on analyzing emission from 
a relatively steady accretion flow, our approach can also 


in the accretion disk (e.g., 


Broderick & Loeb 

2006) or 

jet (e.g., Broderic 

k & Loeb 


20091. Such applications would provide valuable counter- 


parts to potential a strometry of the flari ng region with 
polarimetric VLBI (Johnson et al. 2014) or with near¬ 


infrared int erferometry (e.g., Hamaus et al. 2009 Vincent 
et al.|[20TT). Our method may also be valuable for other 


observations of time-variable structures with sparse visi¬ 
bility data. For example, space-VLBI e xperiments, such 
as RadioAstron ( Kardashe^etaypO^ ) and the planned 
mission Millimetron ( [Wild et al.|2009p will provide many 
baseline pairs with nearly identical length and orienta¬ 
tion from a single space dish to ground stations. 

Moving beyond a simple lagged correlation, rapid rota¬ 
tion of the image could p otentially be used to improve im¬ 
age reconstructions (e.g., Sault et al. 1997). Namely, con¬ 
ventional imaging algorithms assunie that the source is 
static throughout rotation of the Earth and then use the 
Earth’s rotation to increase sampling of the unrotated 
image in the visibility domain (termed Earth-rotation 
synthesis imaging). For the accretion flow of Sgr A*, the 
situation is reversed: the Earth is nearly static for an 
entire rotation of the source. This correspondence sug¬ 
gests that source-rotation synthesis imaging may allow 
an array to achieve better effective visibility coverage on 
timescales of minutes than would be possible with a full 
night observing a static source. Thus, the EHT may be 
capable of rapid snapshot images of an inhomogeneous 
and rapidly-rotating flow, if the rotation curve of the 
source is well-understood. 
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